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We study nuclear high-spin states undergoing the transition to the fully stretched configuration 
with maximum angular momentum I ma x within the space of valence nucleons. To this end, we 
perform a systematic theoretical analysis of non-fully-stretched I ma x — 2 and Imax — 1 fy% seniority 
isomers and 43/2 /"/V intruder states in the A~44 nuclei from the lower- fp shell. We employ 
two theoretical approaches: (i) the density functional theory based on the cranked self-consistent 
Skyrme-Hartree-Fock method, and (ii) the nuclear shell model in the full sdfp configuration space 
allowing for lp-lh cross-shell excitations. We emphasize the importance of restoration of broken 
angular momentum symmetry inherently obscuring the mean-field treatment of high-spin states. 
Overall good agreement with experimental data is obtained. 

PACS numbers: 21.10.Pc, 21.10.Hw, 21.60.Cs, 21.60.Jz, 23.20.Lv, 27.40.+Z 
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I. INTRODUCTION 

The phenomenon of band termination is a splendid 
manifestation of a competition between nuclear single- 
particle and collective motion. At low angular momenta, 
a rotational band is associated with a collective reorien- 
tation of a deformed nucleus in space, with many nu- 
cleons contributing coherently to the total spin. With 
increased rotational velocity, however, the Coriolis inter- 
action causes nucleonic pairs to break, and the decoupled 
nucleons align their individual angular momenta. It of- 
ten happens that breaking relatively few nucleonic pairs 
can give rise to a nuclear state with a fairly large spin. In 
the language of the nuclear shell model, of particular im- 
portance are the "seniority isomers" or "fully-aligned (or 
stretched, or optimal) states," which carry the maximum 
angular momentum I max within the space of valence nu- 
cleons. 

The transition process from collective rotation at high 
spins to a single-particle picture at I ma x is referred to as 
band termination [l| . Terminating bands are common in 
nuclei; they have been observed across the table of the 
nuclides (see recent reviews @, 13)- The nature of the 
termination process strongly depends on the size of the 
valence space. For instance, if only several valence nu- 
cleons are present, the static nuclear deformation cannot 
develop and the collective effects have dynamic charac- 
ter. In this case, I m ax can be generated by breaking 
very few nucleonic pairs, and the transition to the single- 
particle limit is rapid. On the other hand, for deformed 
nuclei having many valence nucleons, the transition to 
the non-collective regime is long and gradual, often in- 
volving many intermediate stages. 

Theoretically, high-spin band terminations and fully 
aligned configurations are often discussed within the nu- 
clear shell model (SM) and/or the self-consistent density 
functional theory (DFT) . The nuclear shell model 0] can 



be applied to nuclei with several valence nucleons out- 
side the magic core. The effective Hamiltonian is exactly 
diagonalized in a subspace of many-body Slater determi- 
nants resulting in correlated wave functions that preserve 
angular momentum, parity, particle number, and - usu- 
ally - isospin. 

For medium-mass and heavy nuclei containing many 
valence nucleons, the tool of choice is the nuclear DFT [f|. 
Here, the nucleus is described in terms of one-body den- 
sities and currents representing distributions of nucleonic 
matter, spins, momentum, and kinetic energy, as well as 
their derivatives. The associated mean fields, obtained 
by means of the self-consistent Hartree-Fock (HF) pro- 
cedure, are usually deformed, i.e., they spontaneously 
break the symmetries of the underlying Hamiltonian. In 
this way, many essential many-body correlations can be 
incorporated into a single product state [1, d, @ . How- 
ever, the price paid for the simple intrinsic picture is 
high: the HF wave function is no longer an eigenstate 
of symmetry operators; hence, the transformation to the 
laboratory reference frame has to be carried out to re- 
store broken symmetries. 

The language in which the nucleus is pictured as a 
wave packet with anisotropic density/current distribu- 
tion is particularly useful for the description of rotational 
motion within the cranking formalism. Here, the many- 
body Hamiltonian H is replaced by a rotating Hamilto- 
nian (Routhian), H u = H — uiJ z , where the rotational 
frequency u> is interpreted as a Lagrange multiplier de- 
termined from the angular momentum constraint. The 
resulting cranked HF method (CHF) can also be used 
to describe the so-called non-collective rotations, i.e., ro- 
tations around symmetry axis 0, 0] • The non-collective 
CHF technique applied in this work has proven to be 
a very efficient and accurate way to generate stretched 
solutions within a self-consistent framework. 

In practice, most applications of nuclear DFT are 
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based on Skyrme energy density functionals optimized 
to various experimental data [E Recently, such an 
approach was applied [EEll to fully-stretched, high-spin 
states associated with the [fy 2 ]i max and [d^J 2 fy J ^']i max 
SM configurations (n denotes the number of valence par- 
ticles outside the 40 Ca core) of 20 < Z < N < 24. 
These studies have demonstrated that those fully-aligned 
states have a fairly simple SM structure, and, therefore, 
they provide an excellent testing ground for both the 
time-odd densities and fields that appear in the mean- 
field description and for SM effective interactions. In 
particular, it was shown that the energy difference be- 
tween the excitation energies of the terminating states, 
E([d~J- 2 ^ l ] Imax ) - E([f? /2 ] Imax ), is a sensitive probe of 
time-odd spin couplings and the strength of the spin-orbit 
term in the Skyrme functional Q, and that with prop- 
erly modifying functionals, the nuclear DFT provides a 
description of the data for stretched states that is of sim- 
ilar quality as the fully correlated SM (Toj . 

The aim of this study is to understand the nature of 
the [/ 7 " /2 ]/ max -i (SecUD and [d^+ 1 ]i max -i (SecEl) 
configurations in the A~44 mass region. These I max — 1 
states, usually referred to as unfavored-signature termi- 
nating states, can be obtained from fully-stretched con- 
figurations by signature-changing particle-hole (p-h) ex- 
citations. Consequently, their structure is sensitive to 
both shell structure and time-odd nuclear fields. We 
demonstrate that non-collective cranking may lead to a 
dramatic violation of rotational symmetry even for the 
cases when the nuclear shape is almost spherical. We 
also analyze in Sec. [V] the I max — 2 weakly collective 
states in terminating normal-parity and intruder struc- 
tures. The associated correlation energy, mainly associ- 
ated with quadrupole effects, is calculated to be large, 
around 2 MeV. 



II. THE MODELS 

The details of DFT and SM frameworks applied in this 
study follow Ref. EU in which the stretched configura- 
tions in the .A~44 mass region have been successfully ex- 
plained. The CHF calculations are carried out using the 
HFODD code of Ref. EH- We employed the SLy4 [l2| 
and SkO [l3[ Skyrme energy density functionals slightly 
modified along the prescription of Refs. @. Without go- 
ing into detail, the modifications concern the coupling 
constants of s 2 and s-As terms giving rise to the time- 
odd spin mean fields. This has been done by constraining 
the functionals to the empirical spin-isospin Landau pa- 
rameters. In addition, the strength of the spin-orbit in- 
teraction has been reduced by 5% from the original SLy4 
and SkO values. As discussed in Ref. jlCj |. the DFT re- 
sults for the T—Q states in N—Z nuclei have to be cor- 
rected for the isospin breaking effects. Here we assume 
that the isospin correction weakly depends on angular 
momentum; hence, it has been neglected when discussing 
energy differences between high-spin states. 



During the last decade, extensive SM calculations for 
the light /p-shell nuclei have been conducted using differ- 
ent SM interactions EE EE USE!. In m ° St CaS6S ' SM 
reproduces well the excitation energies of normal-parity 
and intruder states, as well as transition rates. In general, 
the intruder d^^f™^ 1 states are predicted and observed 
to be more collective than the fy 2 structures, i.e., the as- 
sociated in-band E2 rates are significantly greater within 
intruder bands. The SM results in this work are obtained 
using the code ANTOINE [lj| in the sdfp configuration 
space limited to lp-lh cross-shell excitation from the sd 
shell to the fp shell. In the /p-shell SM space we took 
the FPD6 interaction 12011. The remaining matrix ele- 
ments are those of Ref. [21j . As compared to the earlier 
work [14], the mass scaling of the SM matrix elements 
was done here consistently, thus reducing the sd interac- 
tion channel by ~4%. 



III. UNFAVORED-SIGNATURE TERMINATING 
f? /2 STATES 

We begin with [/ 7 ™ 2 ]; ro „-i configurations in 20 < Z < 
N < 24 nuclei. Figure [1] displays the energy difference 
between the stretched I=I m ax and the lowest I m ax — 
1 states. It is gratifying to see a qualitative, and - in 
most cases - quantitative agreement between SM and 
experimental data. 

Within the naive non-collective cranking, the unfa- 
vored Imax — 1 states can be obtained by either inverting 
the signature of a single proton (ir) or a single neutron 
(i>). The energies of those states, £'j CHF ^' > and Ej CHF ^, 
are displayed in Fig. [T] (top) with respect to the energy 
^(chf) q £ ^ e stretched configuration. It is seen that 

J max 

agreement between CHF and experiment is rather poor. 
In particular, the strong particle number variation in 
the energy difference is not reproduced by CHF. This 
discrepancy has its origin in the spontaneous violation 
of rotational invariance by the mean-field solutions, in 
spite of the fact that the underlying CHF states are al- 
most spherical in all the considered cases. Indeed, since 
the cranking procedure only constrains an expectation 
value of the angular momentum projection, CHF states 
contain components with different angular momentum. 
In the case of unfavored I ma x — 1 states, two compo- 
nents are expected to dominate: spurious reorientation 
mode \I m ax\ Imax — 1) and physical stretched configura- 
tion \Imax 1; Imax !)• 

The appearance of spurious components in partially- 
aligned cranking configurations is well known [22| . A 
classic example is the cranking treatment of two iden- 
tical nucleons in a spherical single- j shell. While the 
stretched configuration with J z =M=2j — 1 can be asso- 
ciated with the stretched state having I max =2j — 1, the 
cranked J z =2j — 2 solution is simply a magnetic sub- 
state of Imax (the Imax — 1 state does not exist in a j 2 
configuration as a result of the Pauli principle). In or- 
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FIG. 1: The energy difference between 

states in /™/ 2 configurations in A~44 nuclei. Dots and cir- 
cles mark experimental data and the SM results, respectively. 
Top: CHF-SLy4 results for unfavored signature terminating 
states corresponding to neutron (filled triangles) and proton 
(open triangles) signature inversion. Bottom: CHF-SLy4 re- 
sults including the angular momentum correction calculated 
using prescription A (filled triangles) and B (open triangles). 



der to remove spurious components, angular momentum 
needs to be restored. Since the number of Af-scheme con- 
figurations around I ma x is very limited, in our cranking 
analysis we resort to an approximate projection scheme. 
In the following, the adopted angular momentum restora- 
tion procedure is discussed for the case of 43 Sc . 

Due to a near-sphericity of the CHF solutions, single- 
particle (sp) states can be labeled using the angular mo- 
mentum projection to. The optimal state \I m ax] Imax) 
can be viewed as a CHF vacuum around which p-h ex- 



citations are built. The \L 
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max 7 ^max 



1) spurious state 
can be obtained by acting on the vacuum with the low- 
ering operator For 43 Sc, which in a SM picture has 
one proton and two neutrons in a / 7 / 2 shell, I max =19/2~~ 
and the (unnormalized) spurious state can be written as 

\Imax>Imax 1) 2\^3 (X^ ^j^v 5/2 \Imax > Imax) ~1~ 

% 5 ^ 2 a 7T r/2\I ma x\Imax) = a\u) +b \ir) , (1) 



TABLE I: Maximum spin states I ma x an d squared unnormal- 
ized expansion amplitudes for fy 2 an d d^^f™^ 1 configura- 
tions in j4~44 nuclei. In the latter case only N^Z nuclei are 
considered. 
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where a\ m (t=tt or v) represents a particle in the /V/2 
shell carrying magnetic quantum number to. By assum- 
ing that the valence coupling scheme ([1]) can be extended 
to the CHF case, the \I m ax] Imax — 1) state can be repre- 
sented as a unique combination of CHF solutions | u) and 
1 7r) corresponding to the lowest neutron and proton p-h 
signature-changing excitations. 

The mixing coefficients a and b introduced in Eq. (fT]) 
for 43 Sc can be calculated in a similar manner for any 
nucleus. They are displayed in Table [T] Assuming that 
contributions from other shells are small, the physical 
state \I ma x — 1; Imax — 1) can be represented by the or- 
thogonal combination: 



Uln 



1) 



(2) 



This 2x2 configuration mixing can be dealt with using 
two slightly different methods. In the first variant (called 
A in the following), one requires that the spurious solu- 
tion is degenerate with respect to the CHF optimal state 
having energy Ej CHF ^ . This assumption leads to a sim- 
ple expression for the energy of the unfavored-signature 
terminating state 
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In the second variant (called B), the mixing coefficients 
are taken directly from Table [J and the corresponding 
energy can be written as 



E 



(CHF,S) 



2p (CHF,jr) 



a J E 



b 2 E (W 



b 2 



(4) 



Note that this method cannot be applied to N—Z nuclei 
where a 2 —b 2 . Moreover, method B does not guarantee 
that the energy of the spurious state is degenerate with 
the optimal state. 



The energy differences between /„ 



and I m ax~l 



states obtained in variants A and B are displayed in the 
lower panel of Fig. [TJ It is seen that the the effect of 
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symmetry restoration is large. In particular, the ener- 
gies corrected for the angular momentum mixing follow 
fairly accurately experiment and SM. Moreover, methods 
A and B give fairly similar results, although method B 
is slightly closer to the data. Similar results were also 
obtained in the CHF-SkO variant. 



IV. UNFAVORED-SIGNATURE TERMINATING 
d 3~/2/7/V STATES 

We now consider the [^/"^'lu.-i intruder config- 
urations. Let us recall that in this case the [d^ 2 f"^ 2 1 ]i max 
configurations are uniquely defined only for N ^ Z nu- 
clei. Indeed, as discussed in Ref. [lfj, the CHF solutions 
for N—Z systems violate isobaric symmetry and can no 
longer serve as reference states. Hence, we shall limit our 
considerations to 20 < Z < N < 24 nuclei. Moreover, in 
the following, we will consider only configurations involv- 
ing lp-lh proton excitation across the Z=20 gap. The 
latter assumption concerns the N— Z=\ nuclei in which 
the neutron cross-shell excitations can also give rise to 
aligned I max — 1 states. It is worth mentioning that the 
lp-lh neutron excitations in N—Z—l nuclei are higher 
than the lowest lp-lh proton excitations. Moreover, they 
do not mix at the level of CHF, leading to a severe isospin 
symmetry violation. Restoration of the isospin symme- 
try for these cases is, however, beyond the scope of this 
work. 

Under these assumptions, within the non-collective 
cranked CHF picture, there are three obvious ways of 
changing the signature quantum number of an intruder 
state. Namely, one can invert the signature of a single 
iV/2 proton ( 1 7r) ) , a single f-j/2 neutron (\v)), or a proton 
d 3 / 2 hole (|7f)). In this case, SM yields two independent 
(Imax — l)i (*=lj2) low-lying solutions. In CHF, both 
these solutions are polluted by the presence of spurious 
\Imax;Imax ~ 1) components. 

Since the CHF solutions for the aligned intruder states 
are nearly spherical, we can employ the same technique 
as proposed earlier for the /™, 2 stretched states to ap- 
proximately restore the rotational invariance. The only 
difference is that in addition to active f-j/2 particles, one 
has to consider active d^/2 proton holes (denoted as 7f in 
the following). As usual, the \I m ax, Imax — 1) is obtained 
by acting with i_ on the \I m ax, Imax) CHF vacuum. For 
the representative case of 42 Ca one obtains: 



Im,n.xilma.x 1) 



\/l2 0,^ g 5 /2 \ Imax ) Imax) 

+ \fl ajj. g/2^?>" 7/2 \Imax i Imax) 
a W l/2 a *3/2|-'maa:) Imax) 

ee a\v) +b\n) +c\tt) . (5) 



Similar calculations can be carried out for all TV ^ Z 
nuclei. The resulting mixing coefficients are collected in 
Table |H 



The intruder case represents a 3 x 3 mixing problem. 
Two different analytical approximate projection tech- 
niques have been developed [23|. In the first method 
(called A), we assume real mixing matrix elements and 
require the eigenvector (a, b, c) of Eq. ([5]) and Table U 
to correspond to zero energy mode relative to the CHF 
energy Ej CHF ^ of the optimal solution. By introducing 



E 



(CHF, a) 



■ Ej CHF \ where a=v (or 1), 7r (or 2), and 



tt (or 3), the energies of physical solutions relative to the 
CHF optimal state i?l CHF ^ are: 

J max 
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In the second method (called B), we admit complex 
mixing amplitudes. In this case we set the spurious mode 
to zero but do not require the corresponding eigenvector 
to be equal (a, 6, c). This procedure leads to exactly the 
same set of formulas ©-(0) as above but with V a p — 
^Je a ep (see Ref. [23( for further details). 

The SM and CHF results corresponding to the lowest 
{Imax — l)i states are shown in Fig. O The agreement 
between SM and experiment is again excellent, except 
for the T—l/2 nuclei 43 Sc and 47 V where theory overes- 
timates the energy difference between I max and I m ax — 1 
intruder states. In order to understand this apparent 
discrepancy, Fig. [2] also shows the SM results for ex- 
cited (Imax — 1)2 states in T—l/2 nuclei. It turns out 
that for 47 V the calculated 33/2^ level has a neutron in- 
truder character, while it is the second 33/2^ state in 
which the proton intruder configuration dominates. For 
45 Ti, the calculated 31/2^" 2 states have a mixed proton- 
neutron intruder character with a slight preference for a 
proton configuration, and the same is true for the 25/2^ 2 
states in 43 Sc. Clearly, the energetics of (Imax — l)i and 
(Imax — 1)2 states in T=l/2 nuclei strongly depends on 
the mixing between proton and neutron intruder states, 
i.e, the isospin dependence of the cross-shell sd <-> fp 
interaction. 

The CHF calculations displayed in Fig. [5] were carried 
out for the SLy4 (top) and SkO (bottom) Skyrme pa- 
rameterizations. It is interesting to see that predictions 
for ^3/3/7/2 intruder states strongly depend on the en- 
ergy density functional. The mean empirical value of 
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44 Ca 44 Sc 45 Ti 



FIG. 2: The energy difference between I ma x and {I ma x — l)i 
states in d^hf?^ 1 configurations in yl~44 nuclei. Dots and 
circles mark experimental data and the SM results, respec- 
tively. The asterisk symbols (SM') show the SM results for 
the (Imax — 1)2 states in T=l/2 Sc, Ti, and V isotopes. The 
CHF results with the SLy4 and SkO functionals are shown in 
the upper and lower panel, respectively. The lowest cranked 
solution (|7f)) is marked by squares while open (filled) trian- 
gles mark the results obtained within the approximate angu- 
lar momentum projection A(B) are marked with open (filled) 
triangles. 



AE = E(I max ) - Ei( I max - 1), averaged over all nu- 
clei considered, is A_Eexp ~ 0.990 keV. For the lowest 
CHF solutions corresponding to fr-type configurations 
(marked by squares in the figures), the calculations yield 
Ai?SLy4 w 0.330 keV and AE S ko » 0.730 keV. The cal- 
culated corrections are of similar size for both Skyrme 
functionals, and in all cases they improve the agreement 
with the data. Within method A, the CHF curves are 
shifted up, on the average, by SAEshyi ~ 0.090 keV 
and (5Ai?sko ~ 0.140 keV. The corrections obtained in 
variant B are SAEshyi ~ 0.210 keV and SAEskO ~ 
0.270 keV. 

In variant B, the average splitting calculated in CHF- 
SkO is fairly close to the experimental value, A_Bgko ~ 
A-Eexp ~ 1 MeV. However, clear discrepancies in the 
isotopic and isotonic dependence are seen. Applying pro- 



jection method B to CHF-SLy4 yields an almost constant 
offset of A~B EX p - A~E S Ly4 ~ 0.450 keV. Indeed, by shift- 
ing the theoretical results up by ~ 0.450 keV, one repro- 
duces surprisingly well the empirical isotopic and isotonic 
dependence, including the cases of 43 Sc and 47 V. It is im- 
portant to point out that both methods A and B are free 
from adjustable parameters. In both cases we have sim- 
plified the picture by limiting the size of the configuration 
space to three states, and in both cases we have forced 
the spurious mode to have zero energy. 



FAVORITE SIGNATURE TERMINATING 
I max - 2 STATES 



In the standard picture of band termination, the 
Imax ~ 2 configurations contain some quadrupole collec- 
tivity. That is, these states have small quadrupole de- 
formations that give rise to a small collective rotational 
component in the wave function. One can therefore ex- 
pect that the energy difference between Imax and Imax— 2 
states should mainly depend on time-even nuclear multi- 
pole fields associated with the nuclear deformability. The 
energy of the last quadrupole transition within the termi- 
nating sequence, AE 2 =E(I max ) - E(I max - 2), is shown 
in Fig.[3]for both fy 2 (top) and d^^f™^ 1 (bottom) struc- 
tures. The agreement between SM and experiment is ex- 
cellent. In 42 Ca, 43 < 44 Sc, and 45 Ti, the B(E2) rates for 
the I m ax—>{Imax — 2)2 transitions are significantly larger 
than those for the yrast I m ax—>(Imax ~ 2)i transitions. 
In these cases, the quadrupole strength is fragmented 
and the band structures cannot be easily identified. It 
is interesting to see that AE2 for the fy 2 configurations 
is significantly smaller (and sometimes even negative) as 
compared to the d^^fj^ 1 values. This is consistent with 
the fact that the quadrupole polarization in the intruder 
states is significantly greater than in the fy 2 structures. 

The mean p-h AJ Z =— 2 excitation energy calculated 
in HF-SLy4 lies ~2 MeV below the data. By performing 
self-consistent cranking, one arrives at a slightly collec- 
tive solution that, on the average, is not too far from ex- 
periment. For the weakly collective, normal-parity struc- 
tures, the CHF theory poorly reproduces the particle- 
number dependence. The situation is significantly im- 
proved in the more collective intruder states where the 
mean-field theory (both in CHF-SLy4 and CHF-SkO 
variants) performs much better. Here, the energy gain 
due to deformation can be as large as ~2 MeV. Of course, 
a big part of the remaining discrepancy between CHF 
results and experiment is due to angular momentum vi- 
olation in the I m ax — 2 states. Unfortunately, for those 
states, the approximate methods of restoring rotational 
invariance discussed in the previous sections do not apply 
because of many states involved (i.e., pronounced config- 
uration mixing). 
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VI. CONCLUSIONS 




SLy4, p-h 



In this work we performed a theoretical analysis of op- 
timal I m ax states and non-fully-stretched I max — 2 and 
— 1 states in the A~44 nuclei from the lower- fp shell. 
Overall, the level of agreement between SM results and 
experimental data for fy 2 seniority isomers and dy\fy + 2 l 
intruder states is excellent. We have shown that CHF so- 
lutions for unfavored-signature terminating states are af- 
fected by dramatic violation of rotational symmetry even 
if the shape of nuclei under consideration is almost spher- 
ical. Approximate methods of restoring rotational invari- 
ance in I max — 1 configurations have been proposed. The 
energy corrections due to angular momentum projection 
can be significant, and their inclusion improve agreement 
with the data. Finally, we investigated the weakly col- 
lective I m ax — 2 members of terminating structures. The 
correlation energy in these states, mainly of a quadrupole 
nature, is fairly large. While for nearly-spherical high- 
spin fy 2 states the CHF method gives only a rough esti- 
mate of the energy splitting between I max — 2 and I max 
states, the agreement is more quantitative for intruder 
configurations. 



42 Ca 42 Sc 44 Sc 44 Ti 46 Ti 47 V 
44 Ca 43 Sc 45 Sc 45 Ti 46 V 



FIG. 3: The transition energy between Imax and {Imax — 2)i 
states in fy 2 (top) and d'^f™^ 1 (bottom) configurations in 
yl~44 nuclei. Dots and circles mark experimental data and 
the SM results, respectively. The intruder (Imax — 2)2 states 
in 42 Ca, 43,44 Sc, and 45 Ti are shown by asterisks. In those nu- 
clei, the B(E2) rates for the Imax^(Imax — 2)2 transitions are 
significantly larger than those for the yrast J ma i->(/ m oi — 2)i 
transitions. The CHF-SLy4 results are marked by a short- 
dashed line. The corresponding average particle-hole excita- 
tion energy is shown by a dotted line. For comparison, the 
CHF-SkO results (dash-dotted line) are given in the bottom 
panel. 
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